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Abstract 

We reexamine the boundary conditions of spin diffusion equations for dirty semiconductor het- 
erostructures with weak Unear Rashba spin-orbit interaction. Doing so, we focus on the influence of 
tangent derivatives of the particle density at the boundary on the magnetization. Such derivatives 
are associated with a spin accumulation in the presence of a density gradient. We show that the 
tangent derivatives enter the boundary conditions and argue that the spin-Hall effect is absent in 
such systems because of this fact. 
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I. INTRODUCTION 



At present there is much interest in spin-charge couphng effects in non-magnetic semicon- 
ductor heterostructures with hnear Rashba spin-orbit interaction, hke the spin accumulation 



effect (see, e.g., [7,8,9,l(f\). While the first two effects can also be observed in homogeneous 
systems the spin-Hall effect manifests itself in a non-equilibrium magnetization at the bound- 
aries of a sample in the presence of electric fields. Thus, investigations of the magnetization 
at the sample boundary require the development of a transport theory, which is also able to 
describe the evolution of particle and spin packets. 

Problems, which are related to the motion of particle packets in semiconductors, have been 
investigated much in the literature (see, e.g., Ref.[ll]]). Drift-diffusion equations are typically 
used to this end. A boundary enters the diffusion equations via boundary conditions, which 
have to be used in seeking solutions to the differential equations. These boundary conditions 
follow often from conservation laws, as it is the case in investigations of charge transport 
processes. There the boundary conditions reduce to the requirement the electric current 
is continuous at the boundary. This condition implies the current vanishes at hard-wall 
boundaries. The structure of the current follows from the conservation of the charge. 

Diffusion equations for the investigation of inhomogeneous magnetizations in semicon- 
ductor heterostructures with Rashba spin-orbit interaction have been derived in a number 
of papers (see, e.g., Refs. ^l2.13^14.15>.16i] ). However, in most cases they have been apphed 
to systems without boundaries. The problem how to formulate boundary conditions for sys- 
tems with hard-wall boundary is a rather controversial one. The conventional conservation 
laws can not be invoked, since the magnetization is not a conserved quantity in the systems 
in question. To circumvent this problem some researchers have taken the point view that 
the spin current can be extracted from the derivative terms in the diffusion equatio n^^i^^ . 
The latter takes the forn>i^ 



in the presence of a constant electric field F [S is the spin density, D is the diffusion constant, 
O is a symmetric tensor of second rank, o is the dyadic product, R = F x Nd^uJs/2 is a 
vector, ujg = AmD/h is a transport coefficient, N = Ne^ is a vector perpendicular to the 
plane of the two-dimensional electron gas, is the derivative with respect to energy, is 
the energy current for spins and Sq = —tN x + tN x Fd^n is the spin accumulation, 
which can be induced either by the constant field i*" or by a gradient of the density n). 
A characteristic feature of this equation is the fact most of its terms depend only on the 
deviation of the magnetization from the spin accumulation So- Thus, the point view that all 
derivative terms are part of the divergence of the spin-current tensor leads to the conclusion 
that the derivative of the magnetization normal to the boundary is proportional to S — So and 
thus zero in the stationary state. Consequently, there is no spin-Hall effect in such systems 
in the presence of a hard-wall boundary. This point of view is supported further by the fact 
that no anomalous Hall-current is induced by injection of particles with perpendicular spin 
in electric fields^. 

However, although diffusion currents have been extracted from diffusion equations for 
centuries it has to be said that this approach is problematic. The diffusion equations contain 
only the divergence of the spin-current tensor, so terms proportional to a curl can not be 




,g5'6 Qj- ^i^Q spin-Hall 



dtS + flo{S-So) + Rx{S- So) - DAS 

+ D{F ■ V)d,S - iu,{N X V) X (5 - 5o) = -dj. 



(1) 
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retrieved unambiguously by this procedure. Moreover, the interpretation of terms as part of 
the divergence of a tensor of second rank is also ambiguous. The sixth term on the left hand 
side (Ihs) of Eq.(ITI), e.g., can also be decomposed into the gradient of a scalar field and the 
curl of a vector field. This raises the question, whether there is another procedure, which 
can be used to determine boundary conditions for diffusion equations. 

Recently, such a method has been developed in Ref.[l3]. Starting from the notion, that 
the boundary conditions for the system in question can be expressed in the form 

en ■ ^Pa = Baf3p(i, (2) 

which connects the normal derivative e„ ■ Vp^ of the densities pa = {n, S) with the densities 
themselves, they calculate the coefficients Ba/s from the integral equations for the densities 
(A summation with respect to double indices must be performed in all equations.). Doing 
so, they use the fact that the integrals are determined by two scales near the boundary, by 
the Fermi wave-length and the mean free path. The first length is the distance from the 
boundary needed to approach a constant particle density and the second the distance needed 
to feel the impact of disorder. Thus, transport is quasi-ballistic in the intermediate region. 
After a tour the force through rather complicated integrals they find that the coefficients 
coupling spin and charge vanish in the intermediate region and that those coupling spin 
and spin are non-zero. This fact has a very important consequence. The derivatives of the 
magnetization do not vanish at the boundary in this case, since the magnetization tends 
to approach the value So in the stationary limit (see Eq.(0)). Therefore, spin-Hall effect is 
obtained. 



The method developed in Ref. [17[ has proven to be very useful in the investigation of the 



s prov 

boundary conditions. The results of Ref.[l7|, however, are neither in line with those of pre- 
vious investigations of the structure of the diffusion equationsi^ nor with the investigations 
of the spin-current tensor. Already a glance on the spin-diffusion equations (P) reveals we 
should expect that the boundary conditions are not determined by the magnitude of the 
magnetization but by the deviation of the magnetization from the spin accumulation Sq. 
The spin accumulation, however, is determined by tangent derivatives. Such derivatives are 
not contained in Eq.Q. Thus, the ansatz Q is simply not sufficiently general enough to 
accommodate the spin accumulation. To allow for an impact of the spin accumulation we 
should replace the ansatz Q by the new ansatz 

e„ ■ Vpa = BapPp + Aa(3 ■ Vp/3 (3) 

and determine the spin-charge coupling coefficients of the matrix Aajj [e-t is the unit tan- 
gent vector). Boundary conditions of the type (jH)) are not unusual for partial differential 
equations. Indeed, also the boundary conditions for Maxwells equations can be formulated 
in terms of normal and tangent derivatives of fields. 

Motivated by our observation we now calculate the spin-charge coupling coefficients of 
the matrix Aa/3 for systems with weak Rashba interaction in the absence of external fields. 
To this end we focus on the situation in a half space and assume that there is a gradient in 
the charge carrier density tangent to the boundary. Moreover, we assume that the normal 
derivative of the charge carrier density vanishes. This is just the situation in a spin-Hall 
experiment. The spin-spin coupling coefficients are of higher order with respect to the 
Rashba interaction and therefore negligible in the systems in question. 
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II. BASIC EQUATIONS 



The Hamilton operator for the two-dimensional electron gas has the form 

H = ^ + N-{(TXp) + V{x). (4) 
2m 

Here p is the momentum operator, cr is the spin-operator and V{x) is a random potential 
with zero average, Gaussian statistics and mean squared deviation 

{V{x)V{x')) = -^6{x-x'). (5) 
Ittut 

The unit vectors in the 2-d plane are e^. and By. is a unit vector perpendicular to the 
plane. To be specific we consider the half-space y > 0. Thus, derivatives with respect to y 
are normal to the boundary and derivatives with respect to x are tangent to the boundary. 
The labels x,y and z are replaced by the indices 1,2,3 in sums. Latin indices run from 1-3, 
Greek indices from 0-3. The index zero is associated with the particle density and ctq is the 
unit matrix. Double indices have to be summed over in all equations. 

The propagation of single particle excitations is described by the retarded and advanced 
Green's functions. To calculate them we use the fact that the system is translation invariant 
in x-direction in average, so we we can use the Fourier transformation. Doing so, we find 
that the retarded (R) Green's functions is given by 

G^{y,y'\k) = -^(e-^«(^+^') - e-'-\y-y'\) (6) 
in the absence of the spin-orbit interaction, where 



\R = -tslkl-k^+'-^. (7) 

Here hkp is the momentum at the Fermi surface and / is the mean free path. The advanced 
(A) Green's function is obtained from Eq.® by hermitian conjugation. To take into account 
the spin-orbit interaction we use the Dyson equation approach of Ref . [l7l| . The spin-orbit 
interaction is considered as the perturbation in this approach and Eq.® as the leading 
approximation. Thus, the approach ignores the impact of the boundary on the scattering 
time. The first correction with respect to the spin-orbit interaction is given by Gf = 
(^iGf + 52Gf , where 

5iGf (y, y'\k) = -a,^^[e-^-l^-^'l(|l/ - y'\ + ^) - e''-(y+y'\y + y' + -^)] (8) 

AriI Ar Ar 



and 

6,G'^{y,y'\k) = aj^{y - y')[e-'-\y-y'\ - e-^«(^^+^')]. (9) 
ARh 



This correction is sufficient to calculate the matrix Bap for systems with weak Rashba in- 
teraction. Therefore, the consideration has been restricted to this contribution in previous 
calculations^'^. The investigation of spin-charge coupling effects, however, requires to go 
beyond this approximation. The leading contributions to these coefficients result from cor- 
rections of second order with respect to the spin-orbit interaction odd in k in the situation 
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discussed here. Therefore, we also need the odd part of the second order correction to the 
Green's function. It is given by 



G^{y,y'\k)U = ,2.4 yy'^.e-'-^'^y'\ (10) 

A nil 



III. THE BOUNDARY CONDITIONS 

To extract the boundary conditions we use the fact that the densities satisfy the system 
of integral equations 

h f 

= J dxitT{aaG^{x,Xi)a(3G'^{xi,x))pf3{xi) (11) 

in the stationary limit. Here the Green's functions are those in position representation. To 
simplify this expression we expand the integrals at the point x and obtain 

Ti f f dk 

Paix) = Pisix)^— J dyi J —tT{a^G^{y,yi\k)afsG^{yi,y\k)) 



Tl P f dk 

+dypi3{x)^^ J dyi J —{yi - y)tT{aaG^{y,yi\k)a(^G^{yi,y\k)) 



+9^Ppi^)^ J dyi J ^tT{a^G\y,y^\k)ap^^kG^{y,,y\k)). (12) 



The first two terms on the right hand side (rhs) of Eq. p2|) are those considered in Ref . jrj 
The third term is the new term associated with the tangent derivative of the fields. To 
calculate the integrals we choose a point x, which satisfies the conditions kpy ':$> 1 and 
y/l <^ 1. Doing so, we find that Eq.(fT^ can be written in the form 

= -DdyS,, = -DdySy - ^ iV^. ^„ (l3) 

= -DdyS, + ^NuJ,{Sy - Soy), (W) 

where D is the diffusion coefficient and Soy is the spin accumulation induced by the density 
gradient. We would like to mention that the terms proportional to the magnetization in the 
Eqs. fllH|l and (fT^ agree with those discussed in Ref. |l7l|. The second term in the bracket of 
Eq.(|14j). however, is a new term, which couples spin and charge. Eq.(|14j) shows, that the 
normal derivative of the magnetization depends only on the deviation of the magnetization 
from the spin accumulation Sn- Thus, there is no spin-Hall effect in systems with hard- wall 
boundary, as detailed in Ref. [16[. 

We would like to mention that the use of the Fourier-transformed Green's functions sim- 
plifies the integrals strongly. If we use the Eqs.Q-CDI) we obtain a number of integrals, 
which contribute to Eq. ()12|) . However, none of them is as complicated as those consid- 
ered in Ref.^3l; where Green's functions in position space were used. The integrals with 
respect to yi reduce to simple exponential integrals, which can be calculated immediately. 
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The remaining integrands reduce either to functions of the type f{k), to functions of the 
type y'^f{k) exp(— ReA/j?/) or to functions of the type y'^Re{f{k) exp(— Ar?/)) (here f{k) are 
complex functions, which do not cause problems). The latter are those, which are strongly 
oscillating, as discussed in Ref . [l7l| . All of the integrals simplify strongly if the chain of sub- 
stitutions k = kpx, x^ — l= esinh(2;), exp(— 2;/2) = w is used (e = 1/kpl). The intermediate 
scale kp^ <^ y can immediately be recognized after these substitutions and a complete 
asymptotic expansion can be obtained with methods as simple as parts integration. We 
would like to mention, that there is also the new parameter ^/ekpy in the integrals. The 
leading contributions, however, do not depend on this parameter. 



IV. ON THE STRUCTURE OF THE DIFFUSION EQUATION 

If we compare the boundary conditions p3|) and (fT^ with those of Ref.(l6l| we note 
they differ from those of Ref. |lq by a factor 1/2. Although this fact does not affect the 
conclusions of Ref. 3] there is the question what the source of this factor is. A glance at 



Eq.(0) reveals, that the factor 1/2 implies, that only 1/2 of the sixth term on the Ihs of 
Eq.(0) is a part of the divergence of the spin-current tensor. Thus, we can write Eq.(^ in 
the form 

9t5 + O o (5 - 5o) + -R X (5 - 5o) + DivJ = (15) 

where 

DivJ = + D{F ■ V)d,S - ]^uOs{N x V) x (5 - 5o), (16) 

and 

R=R- KjsN X V. (17) 

Eq. (jlSp differs from conventional diffusion equations for the magnetization only in that the 
third term on the rhs is proportional to S — Sq. The term x So is absent in conventional 
diffusion equations, since -RH^o in them. Thus, the presence of this term in the diffusion 
equation gives rise to the notion that the impact of the field differs from that of a density 
gradient, as already discussed in Ref.jl^. However, if we take into account the structure of 
the vectors R and R we can also write this equation in the form 

dtS + no {S- So) + Rx S + DivJ = -d,{j, + ^ujs{N x V) x (AT x F)n). (18) 

This fact shows that the impact of the field on the system in question is the same as the 
impact of a density gradient, at least on the level of the diffusion equations. Thus, the 
fact that only 1/2 of the sixth term enters the spin-current tensor restores the symmetry 
between the response to an electric field and the response to a density gradient. Eq. lfTHjl 
also indicates that the second term on the rhs of Eq. ()18j) is a part of the energy relaxation 
current. The diffusion equations for systems with Rashba spin-orbit interaction reduce just 
to the conventional diffusion equations in this case. 

We stress that the conjecture that the second term on the rhs of Eq. p8|) is part of the 
energy current still needs justification. The present paper clarifies the boundary conditions 
in position space. The energy space is a different subject, which we hope to address in a 
future publication. 
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V. CONCLUSIONS 



In the present paper we have reinvestigated the boundary conditions for spin diffusion 
equations with Rashba spin-orbit interaction. Our investigation shows that there are also 
spin-charge couphng terms in the boundary conditions. They couple the tangent derivative 
of the particle density to the normal derivative of th spin-density. There is no spin-Hall effect 
in systems with weak Rashba spin-orbit interaction in the presence of a hard-wall boundary 
due to the spin-charge coupling terms. 
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